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Abstract 

The anomalous behavior of YbInCu4 and similar compounds is modeled 
by the exact solution of the spin one-half Falicov-Kimball model in infinite 
dimensions. The valence- fluctuating transition is related to a metal- insulator 
transition caused by the Falicov-Kimball interaction, and triggered by the 
change in the f-occupancy. 

I. INTRODUCTION 

The intermetallic compounds of the YbInCu4 family exhibit an isostructural transition 
from high-temperature state with trivalent Yb ions in the 4/^^ configuration to the low- 
temperature mixed-valent state with Yb ions fluctuating between 4/^^ and 4/^^ configura- 
tionsB. The transition is particularly abrupt in high-quality stoichiometric YbInCu4 sampleJl 
with a transition temperature equal to T„ = 42 K at ambient pressure; the susceptibility and 
the resistivity drop at Ty by more than one order of magnitude in cooling, while the volume 
expansion is small, AV/V ~ 0.05. The valence change inferred from AV/V by using the 
usual ionic radii of Yb^"*" and Yb^^ is about Anj ^ 0.1, which is consistent with the valence 
measurements by the L///-edge absorptionlll§. The critical temperature depends strongly on 
external pressure, magnetic field, and alloyinglli. A recent review of the experimental data 
is given in Ref.Q and here we just recall the main points which motivate our choice of model. 

The integer- valent phase (T > T^) is characterized by a Curie- Weiss susceptibilityBi 
with very small Curie- Weiss temperature 6 -C T^,. The Curie constant corresponds to the 
free moment of one magnetic f-hole in a J = 7/2 spin-orbit state with /ie// = 4.53/xb. The 
electrical resistance is large and has a small positive slope; it remains almost unchanged in 
magnetic fields up to 30 TO. In some systems, like Ybi_2,.Y^InCu4, the magnetoresistance is 
slightly negative, while in YbInCu4 (or YbIni_^Aga.Cu4 for x = 0.15) it is slightly positive!. 
The Hall constant is large and negative, indicating a small number of carriersQ i. The 
thermoelectric power has a rather small slope which one finds in a semiconductor with a 
nearly symmetric density of statesi. Recent data on the optical conductivity of YbInCu40 
shows the absence of a Drude peak at high temperatures and a pronounced maximum of the 
optical spectral weight at about 1 eV. The high-temperature ESR data for Gd'^+ embedded 
in YbInCu4 resemble those found in integer- valence semi- metallic or insulator hosts0. Thus, 
the high-temperature phase indicates the presence of a well defined local moment but gives 
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no signature of the Kondo effect. The overall behavior of the high-temperature phase is 
closer to that of a semi-metal or paramagnetic small-gap semiconductor than to a Kondo 
metal. 

The mixed-valent phase (T < T„) behaves like a Pauli paramagnet with moderately en- 
hanced susceptibility and specific heat coefficienti. The electrical resistance and the Hall 
constant are one order of magnitude smaller than in the high-temperature phaseBi. The 
thermoelectric poweri has a very large slope typical of a valence fiuctuator with large asym- 
metry in the density of states. The susceptibility, the resistivity and the Hall constant do 
not show any temperature dependence below T^, which is also typical of valence fiuctua- 
tors. The optical conductivity shows a major change with respect to the high-temperature 
shape. The peak around 1 eV is reduced, the Drude peak becomes fully developed, and an 
additional structure in the mid-infrared range appears quite suddenly below T^S. A large 
density of states at the chemical potential /i is indicated by the ESR data as welllii. Thus, 
the transition at seems to be from a paramagnetic semimetal to a valence fiuctuator. 

In contrast to usual valence-fiuctuators, which are quite insensitive to the magnetic field, 
the YbInCu4 family of compounds also exhibit metamagnetic transitions when T < T^. The 
Yb moment is fully restored at a critical field Hc{T), with a Zeeman energy hbHc comparable 
to the thermal energy ksT^. The metamagnetic transition defined by the magnetoresistance 
or the magnetization datJ gives an H-T phase boundary Hc{T) = H^\Jl — (T/T^,)^. The 
zero-temperature field is related to Ty as kBTy/^sH^ = 1-80. 

To account for these features we need a model in which the non-magnetic, valence- 
fluctuating, metallic ground state can be destabilized by increasing temperature or magnetic 
field. Above the transition, we need a paramagnetic semiconductor with an average f- 
occupancy that is not changed much with respect to the ground state. The correct model 
for this system is a periodic Anderson model supplemented with a large Falicov-Kimball 
(FK) interaction term. The temperature or field induced transition suggests that one should 
place the narrow f-level just above the chemical potential fi. The hybridization keeps the 
f-count finite below the transition, while large f-f correlations allow only the fluctuations 
between zero- and one-hole (magnetic) configurations. The low-temperature phase is close 
to the valence fluctuating fixed point and shows no Kondo effect. However, because of the 
Falicov-Kimball term, there is a critical f-occupation at which there is a transition into the 
high-temperature state with a large gap in the d- and f-excitation spectrum. The rif is 
driven to criticality either by temperature or magnetic field. In the high-temperature phase 
the hybridization can be neglected because the f-level width is already large due to thermal 
fluctuations, and quantum fluctuations are irrelevant. Unfortunately, the above model would 
be difficult to solve in a controlled way, and here we consider a simplified model in which 
the hybridization is neglected at all temperatures. This leads to a spin-degenerate Falicov- 
Kimball model which explains the collapse of the non-magnetic metallic phase at or 
He, and gives a good qualitative description of the high-temperature paramagnetic phase. 
However, the deficiency of the simplified model is that it yields a negligible f-count in the 
metallic phase and predicts a large change in the Yb valence at or He. It is clear that we 
can not obtain the valence fluctuating ground state and maintain the average f-occupancy 
below the transition without hybridization-induced quantum fluctuations. In what follows, 
we describe the model, explain the method of solution, and present results for static and 
dynamic correlation functions. 
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II. CALCULATIONS 



The Hamiltonian of the Fahcov-Kimball modelEl consists of two types of electrons: con- 
duction electrons (created or destroyed at site i by dj^ or di^^) and localized electrons (created 
or destroyed at site i by //^ or fi^). The conduction electrons can hop between nearest- 
neighbor sites on the D-dimensional lattice, with a hopping matrix —tij = —t*/2\^; we 
choose a scaling of the hopping matrix that yields a nontrivial limit in infinite-dimensions0. 
The /-electrons have a site energy Ej, and a chemical potential n is employed to conserve 
the total number of electrons n^-f + n^^ + '^/t + "-/j, = ^tot- The Coulomb repulsion Uff 
between two /-electrons is infinite and there is a Coulomb interaction U between the d- 
and /-electrons that occupy the same lattice site. An external magnetic field h couples to 
localized electrons with a Lande g-factor. The resulting Hamiltonian iJllH 

H = ^(-tij - fiSij)dldj^ + Y^{Ef - fi)flU 

+ U diadiafL'fi<T' + Uff fl]fi^fllfi[ 

i,crcr' i,cr 

-f^sh J2 ^C^dladia + gflfia)- (1) 
i,(T 

The model can be solved in the infinite-dimensional limit by using the methods of Brandt- 
Mielschli. We consider the hypercubic lattice with Gaussian density of states p(e) = 
exp[— e^/t*^]/y7ri*, and take t* as the unit of energy (t* = 1). Our calculations are re- 
stricted to the homogeneous phase. 

The local conduction-electron Green's function satisfies Dyson's equation 

G-(^) = / ^feL de, (2) 

where z is a complex variable and S'^ is the local self energy which does not depend on 
momentum0. In infinite dimensions, S'^ is defined by a sum of skeleton diagrams, which 
depend on the local d-propagator G"^ but not on tij. The exact self-energy functional for the 
FK model is obtained by calculating the thermodynamic Green's function0 of an atomic 
system coupled to an external time-dependent field A'^(r) 

Gl^omir) = -^Tr,f (T.e-/^^--c?.(r')4(r)5(A)) , (3) 
where the S-matrix for the A-field is 

g^^^ ^r^^^~jyTjyr'Xiry)dt(r)d^(r')^ (4) 

and -ffatom is obtained from the Hamiltonian (|l|) by removing the hopping and keeping just a 
single lattice site. The exact solution for ^^[{Am}] at Matsubara frequency iujn = i7iT{2n+l) 
is given by, 

where Wq and Wi are the f-occupation numbers {wi = 1 — Wq, Wq = Zq/Z) andEl 
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Zo(A, ,) = 2e^^l^ n 7-^2e^-/^ R 7"^, (6) 

with 

Z{\, fi) = Zo{\, /i) + 2e-^(^^-'')Zo(A, fi-U). (7) 
The bare Green's function satisfies 

^on = ■ 7^. (8) 

ZUn + /i — A„ 

with A„ the Fourier transform of the external time-dependent field. 

The self-energy functional can now be obtainedEl by using the Dyson equation 

for the atomic propagator, 

K = [G^on]-' - [GT\ (9) 

and eliminating Gonfl-^m}] from Eqs. (|) and @. The mapping onto the lattice is achieved 
by adjusting Gq„ such a way that G^[{Am}] satisfies the lattice Dyson equation (^. 

The numerical implementation of the above procedure is as follows: We start with an 
initial guess for the self energy S°" and calculate the local propagator in (|^). Using (^) we 
calculate the bare atomic propagator Gq^ and find Zq and Z. Next we obtain wq, Wi and 
find from (^). Using Gq^ and G'^, we compute the atomic self energy and iterate to the 
fixed point. 

The iterations on the imaginary axis give static properties, like n/, the f- magnetization 
m/(/i, t), and the static spin and charge susceptibilities. Having found the f-electron filling 
wi at each temperature, we iterate Eqs. (||) to (||) on the real axis and obtain the retarded 
dynamical properties, like the spectral function, the resistivity, the magnetoresistance, and 
the optical conductivity. At the fixed point, the spectral properties of the atom perturbed 
by A-field coincide with the local spectral properties of the lattice. 



III. RESULTS AND DISCUSSION 

We studied the model for a total electron filling of 1.5 and for several values of Ef and 
U. The main results can be summarized in the following way. 
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FIG. 1. Number of the /-holes plotted versus T/t* for U/t* = 4. The Ej/t* increases from 
top to bottom, and is given by -0.7, -0.6, -0.5, and -0.2, respectively. 

The occupancy of the f-holes at high temperatures is large and there is a huge magnetic 
degeneracy. The f-holes are energetically unfavorable but are maintained because of their 
large magnetic entropy. In Fig.(|l]) we show n/ as a function of temperature, plotted for 
U = At*, and Ef/t* from -0.2 to -0.7. Below a certain temperature, which depends on U 
and Ef, there is a rapid transition to the low-temperature phase. The transition becomes 
sharper and is pushed to lower temperatures as Ef decreases at constant U. However, we 
restrict ourselves to continuous crossovers here, since the region with first-order transitions 
leads to numerical instabilities. 

The uniform f-spin susceptibility is obtained by calculating the spin-spin correlation 
functionlll'0 and is given by x{T) = Cnf{T)/T, where C = g\^\J{J + l)/3fcB is the Curie 
constant. 
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FIG. 2. Uniform static magnetic susceptibility of the /-holes plotted versus T/t* for U /t* = 4. 
The values of Ej/t* are the same as in Fig(|l]). The corresponding values of T^/t* are estimated 
from the maximum of given by 0.03, 0.08, 0.15, 0.35, respectively. The Ty increases 

from top to bottom. 

The x{T)/C is shown in Fig. (Q) for U/t* = 4 and for Ef as quoted in Fig.(|l|). The 
Ty is obtained from the maximum of the x{T)/C and the values corresponding to various 
parameters used in this paper are quoted in the caption of Fig. The high-temperature 
susceptibility follows an approximate Curie- Weiss law, but the Curie- Weiss parameters de- 
pend on the fitting interval. 

The interacting density of states pdiyj) for the conduction electrons is shown in Fig.(^ 
for U/t* =4 and Ej/t* = —0.5, and for several temperatures. (The energy is measured 
with respect to /i.) The high-temperature DOS has a gap of the order of U, and the 
chemical potential is located within the gap. Below the transition is small, the correlation 
effects are reduced, and pd{uj) assumes a nearly non-interacting shape, with large pd{p) and 
halfwidth ~ t*. 

The transport properties of the high-T phase are dominated by the presence of the gap, 
which leads to a small dc conductivity with a weak temperature dependence. The transport 
properties of the paramagnetic phase are unrelated to the spin-disorder Kondo scattering 
(there is no spin-spin scattering in the FK model). Below the transition the conductivity 
increases and assumes large metallic values. 
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FIG. 3. Interacting density of states plotted versus Lo/t* for U/t* = 4, Ef/t* = —0.5 
{Ty/t* = 0.14), and for various temperatures, as indicated in the figure. 

The intraband optical conductivity cioj) is plotted in Fig.(^ as a function of frequency, 
for several temperatures. Above T^,, we observe a reduced Drude peak around = and 
a pronounced high-frequency peak around ~ f/. The shape of criuj) changes completely 
across Ty. Below Ty the Drude peak is fully developed and there is no high-energy (intraband) 
structure. However, if the renormalized f-level is close to /i, the interband d-f transition 
could lead to an additional mid-infrared peak. The ratio of the high-frequency peak in 
Fig(^) and the corresponding value of Ty = O.lSt*, is U/Ty = 26. For the same value of U 
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and Ef = -0.7t* (T^ = 0.03t*) we obtain U/T^ = 130, while for Ej = -0.75t* (T^ = 0.02t*) 
we find U/T^ ~ 200 (not shown). 
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FIG. 4. Optical conductivity plotted versus uj/t* for various temperatures. The U, Ef, and T^, 
are the same as in Fig.(^. 

If we estimate the f-d correlation in YbInCu4 from the 8000 cm~^ peak in the optical 
conductivity data0, we obtain the experimental value f/ ~ 1 eV. Together with = 42 
kI this gives the ratio f//T„ ~ 200. If we take U/t* = 4 and adjust Ef/t* so as to bring 
the theoretical value of in agreement with the the thermodynamic and transport data on 
YbInCu4, we get a high-frequency peak in cr^u) at about 8000 cm~^, 6000 cm~^, and 1500 
cm-\ for Ef = -0.75t*, Ef = -0.7t*, and Ef = -0.5t*, respectively. 




FIG. 5. The f-electron magnetization mj is plotted as a function of h/t* for various tempera- 
tures. The U, Ef, and Ty, are the same as in Fig.(^). 

The f-electron magnetization rnf{h) is plotted in Fig.(|^) versus reduced magnetic field 
h/t*, for several temperatures. Above the characteristic temperature T* ~ Ty/2, the mf{h) 
curves exhibit typical local moment behavior. Below T* we find a metamagnetic transition 
at a critical field He, the mf{h) is negligibly small below He and the local moment is fully 
restored above Hr. 
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FIG. 6. Normalized critical field is plotted as a function of reduced temperature T/T* for 
several values of Ef/t* and U/t*. The full line represents ^Jl — {T /T*)'^ and T* = 

Taking the inflection point of the mf{h) curves, calculated for several values of U and Ef, 
as an estimate of Hc{T) we obtain the phase boundary which is shown in Fig.(|^), together 

with the expression H^{T)/H^ = ^1 - {T/T*^. Note, the T* values in Fig.(|) differ by 
more than an order of magnitude, while the ratio ksT* / fisH^ is only weakly parameter 
dependent. 
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FIG. 7. Field-dependent resistivity plotted versus h/t*. The different symbols correspond to 
different temperatures, as indicated in the figure. The U and Ef are the same as in Fig.(^. 

The metamagnetic transition is also seen in the field-dependent electrical resistance 
R{h,T) which is plotted in Fig.(|^) as a function of h/t*, for several temperatures. A sub- 
stantial change in the R{h, T) across T* or He is clearly seen. 
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IV. SUMMARY 



From the preceding discussion it is clear that Fahcov-Kimball model captures the main 
features of the experimental data for YbInCu4 and similar compounds. The temperature- 
and field-induced anomalies are related to a metal-insulator transition, which is caused by 
large FK interaction and triggered by the temperature- or the field-induced change in the 
f-occupancy. At high temperatures, we find a large gap in Pdi^^)', we expect a similar gap 
in the f-electron spectrum as well. At low temperatures, both gaps are closed, and the 
renormalized f-level renormalizes down to the chemical potential. 

Our calculations describe doped Yb systems with broad transitions but appear to be less 
successful for those compounds which show a first-order transition. The numerical curves 
can be made sharper (by adjusting the parameters) but they only become discontinuous 
in a narrow parameter range. The main difficulty with the FK model is that it predicts 
a substantial change in the f-occupancy across the transition and associates the loss of 
moment with the loss of f-holes. But in the real materials the loss of moment seems to be 
due to the valence fluctuations, rather than to the reduction of Uf. The description of the 
valence fluctuating ground state would require the hybridization and is beyond the scope of 
this work. The actual situation pertaining to Yb ions in the mixed-valence state might be 
quite complicated, since one would have to consider an extremely asymmetric limit of the 
Anderson model, in which the ground state is not Kondo-like, there is no Kondo : 
and there is no single universal energy scale which is relevant at all temperatures 

We speculate that the periodic Anderson model with a large FK term will exhibit the 
same behavior as the FK model at high temperatures. Indeed, if the conduction band and 
the f-level are gapped, and the width of the f-level is large, then the effect of the hybridization 
can be accounted for by renormalizing the parameters of the FK model. On the other hand, 
if the low-temperature state of the full model is close to the valence-fluctuating fixed point 
with the conduction band and hybridized f-level close to the Fermi level, then the likely 
effect of the FK correlation is to renormalize the parameters of the Anderson model. 
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